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ABSTRACT 

About half of the binary systems are close enough to each other for mass to be exchanged between them at some point in their 
evolution, yet the accretion mechanism in wind accreting binaries is not well understood. We study the dynamical effects of 
gravitational focusing by a binary companion on winds from late-type stars. In particular, we investigate the mass transfer and 
formation of accretion disks around the secondary in detached systems consisting of an AGB mass-losing star and an accreting 
companion. The presence of mass outflows is studied as a function of mass loss rate, wind temperature and binary orbital 
parameters. A 2-dimensional hydrodynamical model is used to study the stability of mass transfer in wind accreting symbiotic 
binary systems. In our simulations we use an adiabatic equation of state and a modified version of the isothermal approximation, 
where the temperature depends on the distance from the mass losing star and its companion. The code uses a block-structured 
adaptive mesh refinement method that allows us to have high resolution at the position of the secondary and resolve the formation 
of bow shocks and accretion disks. We explore the accretion flow between the components and formation of accretion disks 
for a range of orbital separations and wind parameters. Our results show the formation of stream flow between the stars and 
accretion disks of various sizes for certain orbital configurations. For a typical slow and massive wind from an AGB star the 
flow pattern is similar to a Roche lobe overflow with accretion rates of 10% of the mass loss from the primary. Stable disks 
with exponentially decreasing density profiles and masses of the order 10 solar masses are formed when wind acceleration 
occurs at several stellar radii. The disks are geometrically thin with eccentric streamlines and close to Keplerian velocity profiles. 
The formation of tidal streams and accretion disks is found to be weakly dependent on the mass loss from the AGB star Our 
simulations of gravitationally focused wind accretion in symbiotic binaries show the formation of stream flows and enhanced 
accretion rates onto the compact component. We conclude that mass transfer through a focused wind is an important mechanism 
in wind accreting interacting binaries and can have significant impact on the evolution of the binary itself and the individual 
components. 

Subject headings: Accretion, accretion disks - Binaries: symbiotic - Hydrodynamics - Methods; numerical - 
Stars: mass loss - Stars: winds, outflows 



1. INTRODUCTION 

Accretion is a very important energy source for a wide va- 
riety of astronomical objects including many interacting bina- 
ries. The mass transfer in interacting binaries usually occurs 
via tidal interaction and Roche lobe overflow (RLOF). How- 
ever, a wide range of interacting binaries belong to a class of 
accreting systems where the primary does not fill its Roche 
surface, and the mass transfer is believed to occur mostly by 
wind accretion. These include many symbiotic binaries and 
even more exotic systems such as massive X-ray binaries and 
microquasars. Wind accretion is more complicated than tidal 
interaction, and the accretion processes in wind accreting bi- 
naries are poorly understood. In this paper, we study gravita- 
tionally focused wind accretion in symbiotic systems with an 
AGB star and a compact accretor where the photosphere of 
the primary does not fill its Roche lobe. 

Symbiotic binaries are important astrophysical laboratories 
for studies of wind accretion because of the wide separa- 
tion of the components, and the ability to study the individ- 
ual components and the accretion processes at many wave- 



lengths ranging from X-ray to radio ( |Karovska et al. 1997 
[2005 , Matthews & Karovska 2006). A typical symbiotic con- 
sists of a mass-losing AGB or a red giant star and a hot accret- 
ing companion, often a white dwarf (WD). The components 
in these systems are assumed to be detached (at least during 
most of the orbital motion) and the compact companion ac- 
cretes mass from the massive wind of the cool evolved star 
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Mass loss is known to be a key factor in the late stages of 
the evolution of an AGB star Matter escapes easily because 
of the low surface gravity. As the intermediate mass stars en- 
ter the AGB phase their atmospheres rise to sizes of the order 
of few AU. These stars become unstable to large amplitude ra- 
dial pulsations and drive strong shock waves that produce an 
extended envelope and a copious stellar wind driven by radi- 
ation pressure on dust grains that form at several stellar radii 
above the photosphere ( |Bowen|fl988[ ). The powerful stellar 
winds also affect the evolution of the ISM by influencing its 
dynamics and spreading elements generated in stellar interi- 
ors. It is believed that interstellar dust has its origin in cool at- 
mospheres and dense winds blowing from these evolved stars. 

Due to massive stellar winds, intermediate mass stars can 
avoid exploding as a supernova and become planetary nebu- 
lae (PN) and p roto-planetary nebul ae (PPN) at the end of the 
AGB (see e.g. Gawryszczak et alT| p002). Mass-losing giant 
stars and PPNs often present asymmetric envelopes and bipo- 
lar outflows that can be explained by wind gravitational focus- 
ing and wind collision dynamics. The presence of bipolar and 
collimated outflows in PPNs has been explained due to a colli- 
sion of stellar winds (Mellema & Frank 1997 , Gonzale z et al.| 
|2004| l and magnetically driven jets (Blac kman et al.|2001 1. 

In the following we show that although the photosphere of 
the AGB star in symbiotic systems does not fill its Roche lobe, 
its extended atmosphere can be significantly disturbed by the 
accreting companion. Mass accretion can be increased due to 
the gravitational attraction from the accretor and the centrifu- 
gal force due to the rotation of the system. This is important, 
because the binary parameters can evolve which will affect 
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the stellar mass loss and wind parameters. Mass and angular 
momentum can be lost from the system or accreted onto the 
companion. The angular momentum loss can result in orbital 
shrinking of the system which would enhance mass transfer 
and can lead to a Type la supernova (Hachisu et al. 1999{ |Ja-| 
lhanara et ai;]p005) |PoJsiadlowski & Mohamed 20071 . The 
evolution of AGB stars in binary systems gives rise to highly 
aspherical winds and can be significantly different from that 
of single stars. 

We use 2-dimensional hydrodynamical time-dependent 
models to investigate the mass transfer driven by a massive 
AGB wind that is captured by a detached companion star. 
We calculate accretion rates assuming wind acceleration at 
different radii of the extended atmosphere (dust formation 
radii) and also consider the case when the extended envelope 
reaches the Lagrangian point Li and the gas flows in a nar- 
row stream towards the secondary. In Section |2] we discuss 
Bondy-Hoyle wind accretion in a symbiotic systems with an 
AGB mass-losing star. We describe our numerical setups in 
Section |3] In Section |4] we present the results of the numer- 
ical simulations and we compare the theoretical models with 
recent observations of symbiotic binaries to improve our un- 
derstanding of the properties of these systems. Finally, we 
discuss the results and future work in Section|5] 

2. BONDI-HOYLE ACCRETION 

Physical processes involved in winds from extended atmo- 
spheres are very complex and only recently a detailed theo- 
retical modeHnghasteen^^ e.g. Bowen & Willson 



T99Tl|Arndt et al.'1997','Mastrodemos & MoiTis .l998„Sandin| 



& Hofner|2004[ |Freytag & Hofner 2008). In the classical the- 
ory of plane-parallel stellar atmospheres the physical structure 
of the atmosphere is determined by the effective temperature 
and surface gravity for a given chemical composition of the 
atmosphere. In the case of the AGB stars the situation be- 
comes more complex. Pulsations in the atmosphere, phenom- 
ena generated in the stellar interior such as convection, and 
especially dust formation influence the mass loss and wind 
parameters. 

For high velocity winds where the pressure can be ne- 
glected, the Bondi-Hoyle accretion rate is given by (see Edg^ 
2004 for a recent review of Bondi-Hoyle-Lyttleton accretion) 



where pb is density at the secondary's position, Vi 



(1) 



v4 + v'g is the wind velocity relative to the accreting object, 

v„ is the wind velocity at the position of the secondary, vb 
is the orbital velocity of the companion, and /?a is the Bondi- 
Hoyle accretion radius given by 

2GMb 



When the mass accretion 
dimensional case we obtain 



(2) 

rate is normalized to the 2- 



Mbh = 2R.,p^v,. (3) 
The mass loss from the unperturbed wind of the AGB com- 
ponent in the orbital plane at the secondary position is 

M^ind = 27TdpBV„. (4) 

Therefore the Bondi-Hoyle accretion ratio, defined as the 
mass accretion divided by the mass loss rate from the primary, 
in our 2-dimensional model is given by 

Mbh 



/bh - 
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(5) 
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Fig. 1 . — Diagram of the accretion geometry around the secondary marked 
by the gray circle. The flow is deflected at the shock wave and most of the 
accretion takes place along the accretion line behind the accretor 

with q denoting the mass ratio, d the binary separation and we 
have used the definition of the orbital Keplerian frequency 



Ob = 



G(Ma + Mb) 
d^ 



1/2 



(6) 



As we will see in Section |4!2) the fraction of the AGB wind 
accreted by the secondary is several times higher than the 
Bondi-Hoyle rate in our simulations due to gravitational wind 
focusing. However, the Bondi-Hoyle theory cannot be di- 
rectly applied to the wind accretion scenario in the case of 
slow winds where the sound speed is not negligible. 

If the gas pressure is included in the Bondi-Hoyle scenario 
( |Bondi 'l952 ; Theuns & David| 1992) 1, the mass accretion ratio 
in the orbital plane becomes 



2j2 



/b: 



n^d 
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where ^ is the Mach number 



(8) 



and Csis the sound speed. The accretion radius is defined as 

2GMb 



(9) 



A diagram of the flow structure in the orbital plane is shown 
in Fig.[T| The stellar wind flows from the left and impacts the 
shock wave and is deflected towards the accreting companion. 
In the Bondi-Hoyle picture, the gas reaches the accretion line 
where the transversal velocity component is neglected and the 
gravitationally bounded material is accreted. In our wind ac- 
cretion model the AGB photosphere does not fill its Roche 
lobe and the mass loss is carried out in the form of a slow 
wind from the region where dust grains are condensed. 

The wind in our models is characteri zed by a typical A GB 
mass-loss rate M ~ 10 ''-10 Mq yr~' (|Ryde et al.|2000|) and 
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the unperturbed initial velocity profile is given by the isother- 
mal wind solution. In our slower wind models the Bondi 
accretion radius is comparable to the stellar separation and 
Bondi-Hoyle theory is not valid. Nevertheless, we will use the 
Bondi-Hoyle rates to compare with the mass accretion rates 
resulting from our simulations. 

3. NUMERICAL MODEL 

Two different problems in the context of mass loss in late 
stages of stellar evolution are investigated. In the first case the 
primary's extended envelope fills its Roche lobe and the gas 
flows in a narrow stream towards the secondary. This could 
be the case when the companion is in periastron. We also use 
this model to test our numerical model. Mass transfer in the 
second model is driven by a massive wind that is captured by 
a detached companion star. This is an important mechanism 
for the formation of accretion disks in detached interacting 
binaries, symbiotic stars and massive X-ray binaries (see e.g. 
ILivio & Wamer|1984 , Soker 2004 ). 

We follow the gas flow in the orbital plane of the system us- 
ing 2-dimensional hydrodynamical simulations that allows us 
to resolve the large density contrasts close to the secondary. 
We limit our wind accreting model to the orbital plane of the 
binary system where the wind is focused forming an equato- 
rial outflow that has a significant density enhancement com- 
pared with the regions further polewards ( Walder et al. 2008| l. 
In Fig. |2] we show a diagram of the Roche potential in the 
equatorial plane for a mass ratio q - 1. 

The basic equations of hydrodynamics describing the evo- 
lution of the density and velocity field are: 



J+V-(pv) = 
ot 

d\ 1 

— + (v ■ V)v = — Vp - Vd), 

at p 



(10) 
(11) 



where p is the density in the orbital plane, v the velocity of the 
fluid, p the pressure and (D the gravitational potential. Both 
stellar gravity potentials were approximated as point masses. 
To avoid numerical problems close to the stars the potential is 
softened using the function: 



Vlr-rAp + 62 



Vlr-rBp + e2 



(12) 



where the smoothing length e was set to 0.01 length units in 
our simulations with standard resolution. The unit of length 
in our calculations corresponds to the separation between the 
stars. In our higher resolution simulations the number of grid 
zones in the accretion region was ~ 200. The surface of the 
compact accretor was not resolved in our models. 

3.L Numerical code 

Our code is based on FLASH ( [Fryxell et al.|200()l l, which is 
a fully parallel block-structured Adaptive Mesh Refinement 
(AMR) implementation of the Piecewise Parabolic Method 
(PPM) algorithm (Woodw ard & Col ella| [T984t |Colella~& 
Woodward 1984 1 in its original Eulerian form ^. ITie code has 



been extensively tested in various compressible flow problems 
with astrophysical app lications (see e.g. |Calder et al.|2002 de 
IVal-Borro et al.|2006| l. 

For numerical convenience we introduce dimensionless 
units, where the binary separation d is taken as the unit of 

^ The FLASH code is available at |http://www.flash.uchicago.edu/| 



length. The unit of time is calculated from the orbital angular 
frequency Qb of the system. The orbital period of the system 
is then defined by 

Pb^2K. (13) 

Each simulation was run for about 10 orbital periods, when 
the system has reached a quasi-equilibrium state. The evolu- 
tionary time of the results discussed in Section[4]will be given 
in units of Pb- 

We normalize all physical quantities for numerical conve- 
nience. The gravitational constant G is set to one as is done 
frequently in relativistic calculations. Eq.|6]can be considered 
as a normalization condition for the unit of mass to be the sum 
of the stellar masses M - Ma + Mb- 

Our implementation of the FLASH code uses polar coor- 
dinates and was run in both the coordinate system corotating 
with the angular speed of the system and in the inertial frame. 
In addition to the usual hydrodynamic quantities, our model 
includes the gravity of both stars, and the Coriolis and cen- 
trifugal forces. The code is based on release 2.5 of FLASH 
with customized modules for the equation of state and gravity 
forces that explicitly conserve angular momentum transport 
( |Kley| |l998 ). This is particularly important when large den- 
sity gradients are present in the wind. The Coriolis forces 
were treated conservatively as described by |Kley| ( |1998 b. A 
Courant number of 0.7 was used in the simulations. We use 
an equation of state for an ideal gas with ratio of specific heats 
between y = 1 - 5/3 to represent more realistic cases that in- 
clude radiative cooling, and a modified isothermal Riemann 
solver ported from the AMRA code ( Plewa & Muller 2001) . 

We ignore magnetic fields, radiation transfer or explicit 
Navier-Stokes viscosity in our models. FLASH has a small 



numerical viscosity that has been estimated by de Val-Borro 
et al.| ([2007). Non-axisymmetric effects on the tidally en- 



hanced wind may be caused by magnetic fields. It is likely 
that stellar winds in AGB stars are primarily radial since the 
extended envelope cannot have significant rotation (Soker & 
Zoabi 2002 ). The MHD collimation of the wind is a very effi- 
cient mechanism and can be more important than gravitational 
focusing (Garcia-Segura & Lopez 2000; Garcia-Segura et al. 
|2001| l. However, MHD effects are not included in our calcula- 
tions since there are presently no constraints on the magnetic 
fields in symbiotic systems. At the AGB phase, it is unlikely 
that a globally ordered magnetic field strongly affects the dy- 
namics of the wind. 

3.2. Mass accretion onto the secondary 

The formation of an accretion disk can be affected by accre- 
tion onto the secondary. In the case of a reflecting boundary 
the stream bounces off the stellar surface and orbits the ac- 
cretor However, if no mass is allowed to accrete onto the 
secondary the orbiting disk continues to increase as the sim- 
ulation continues. We use a mechanism that allows us to re- 
move gas from the vicinity of the secondary. The accretion 
is accounted for by removing some mass from the region de- 
fined by |r - TbI < race and adding it to the mass of the star 
to calculate the gravitational forces. The goal was to obtain a 
quasi-stable configuration and to be able to estimate the accre- 
tion rate onto the accretor The size of the accretion region race 
is defined as a fraction of the Hill radius of the mass accretor. 
In most of our simulations the accretion radius is race - 0.1/?h, 
where Rn is the Hill radius of the secondary 



Mb 
3Ma 



1/3 



(14) 
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The mass is rem oved from the disk after each time step us- 
ing the expression ( Giinther & Kley|2002| i 

Ap = min (0, /Af) max (0,p - p^y) , (15) 

where / is a constant fraction of the order unity, Af denotes the 
timestep and pav is th e average density in the region r^cc < |r- 
TbI < 2racc following Peplinski et al. (2008 1. The momentum 
of the accreted material is removed from the system. 

We calculate the mass accretion ratio as the ratio of mass 
captured by the companion to the mass loss rate from the pri- 
mary averaged over 100 timesteps 



^ ^ / Mace \ 



(16) 



[Theuns et al.| ( |1996l l have found accretion efficiencies around 
2% for binaries with mass ratio q -2. 

3.3. Roche lobe overflow 

Some of the most dramatic phenomena in binary systems 
occur when an AGB star fills its Roche lobe. In this case, 
a dynamically unstable mass transfer to the companion takes 
place that can lead to the formation of a common envelope and 
a very fast evolution of the binary system. We model the in- 
teraction of a tidal Li stream with an accretion disk around 
a compact star using an adiabatic equation of state and an 
isothermal solver to compare the results with the wind accre- 
tion simulations and test our numerical code. 

The density of the stream flowing from the Lagrangian in- 
ner point has the form. 



P =Po + 



Ps -PO 



cosh((0/<^o)'") 



(17) 



where po is the initial unperturbed density, ps is the density of 
the stream, is the azimuthal angle of the inner Lagrangian 
point and m is a coefficient of order 20 that gives the steepness 
of the stream. 

The initial radial velocity was of the order 1 % of the orbital 
velocity of the system and directed towards the mass accreting 
companion. We have checked that diff'erent initial velocities 
at the Li point do not affect om- results. According to the an- 
alytical calculations of |Lubow & Sh u ( 1975 ), the tidal stream 
leaves the first Lagrangian point at an incidence angle of ap- 
proximately 28° with respect to the line joining the centers 
of the two stars. Since the unit of density po drops out of 
the equations of motion, we may normalize it to any specified 
density. It is useful to set the stream density to a fraction of 
the total mass of the system to be able to compare accretion 
rates to the initial mass. 

We used two different boundary conditions at the inner edge 
of the grid. The stream flow and formation of an accretion 
disk depend on the interaction of the tidal stream with the 
surface of the accreting star. Therefore it is important to test 
different boundary conditions at the surface of this star. We 
used reflecting boundary and allowed accretion onto the star 
using the prescription defined in Section 3.2 and an absorbing 
boundary where the radial gradients of the velocity pressure, 
and density are zero at the inner boundary. Simulations us- 
ing an absorbing boundary allow the stream to flow into the 
accretor without generating reflected waves. We assume that 
the computational domain is filled with tenuous gas with the 
temperature profile of the wind. Initially, the density of the 
ambient material is 10 orders of magnitude smaller than the 
density at the surface of the mass-losing star. 




Fig. 2. — Roche potential on the orbital plane showing the location of the 
Lagrangian points for a mass ratio g = I. 




Fig. 3. — Grid geometiy of the RLOF simulations in polar coordinates. The 
computational domain is centered on the accretor and a stream introduced 
over a few cells at the inner Lagrangian points. 



The hydrodynamic flow was computed in cylindrical coor- 
dinates centered on the accreting companion and rotating in 
the reference frame of the binary system. Most of the sim- 
ulations were run on a grid with 128 radial zones times 384 
angular zones and 3 additional levels of refinement. The radial 
spacing in our base grid is uniform. Fig. [3] shows a diagram 
of the grid geometry and boundaries in the RLOF setup. The 
computational grid was centered on the accretor and a stream 
is initially introduced over a few cells at the inner Lagrangian 
point. Some tests were run with twice the linear resolution to 



check the convergence of the results (see Sect. 4.2.4 1. 

Two different equations of state were used in our model. 
One version uses an isothermal Riemann solver and the adia- 
batic models use an ideal gas equation of state with polytropic 
index y in the range 1 to 5/3. The isothermal code was used 
to explore the parameter space because of its relative speed 
without the need to calculate the cooHng and heating terms 
and evolving the energy equation. 
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Fig. 4. — Radial velocity profile of the unperturbed wind for a temperature 
of lO^K. The dotted line corresponds to the escape velocity from a 1.2Mo 
star. 

3.4. Wind accretion 

In this section we describe our numerical setup to study 
wind accretion in detached symbiotic binaries where the AGB 
component does not fill its Roche lobe. Many detached inter- 
acting binaries, including symbiotic binaries and X-ray bina- 
ries show evidence for much higher mass transfer rate than 
predicted by models of spherically symmetric wind. The en- 
hanced accretion rates onto the secondary can be explained by 
wind focused towards the compact companion. 

3 .4. 1 . Wind model 

The initial spherically symmetric wind structure is obtained 
by solving the momentum equation 



dv, 1 dP 

Vr 1 

dr p dr 



GMa 



(18) 



where v,. is the radial component of the wind velocity and 
Ma the mass of the central star The solution for an isother- 
mal wind is shown in Fig. [4] As the wind propagates across 
the grid it is accelerated by thermal pressure and its velocity 
increases. It is possible that AGB winds are not completely 
spherically symmetric due to i nternal processes in the star or 
the influence of the secondary ( Dorfi & Hoefner|1996 1. 

In Fig. |5]we plot the temperature profile used in our wind 
model. During the AGB phase it is likely that the winds are 
primarily radial since the massive stellar atmosphere cannot 
have significant rotation. Our temperature distribution has a 
radial profile oc r and is similar to the one given by the 



AGB models by |Ivezic & Elitzur| ( |1997) l. Fig|5] shows the 
temperature distribution after 6 orbits in our locally isother- 
mal simulations. 

The binary systems considered here consist of a 1.2 M© 
AGB star and a 0.6 Mq companion. The orbital separation d 
ranges between 10 AU and 70 AU. The initial temperature of 
the wind is between 3000 K and lOOOOK and the velocity of 
the wind on the surface of the AGB star is 2 - 15km s ' . Both 
temperature and the radial velocity of the wind are constant at 
the dust formation radius for a given wind model. 

The computational domain was centered on the primary and 
the secondary was fixed on a circular orbit neglecting the hy- 
drodynamic effects. Our grid extended out to twice the bi- 
nary separation in our models which is enough to avoid mass 



Fig. 5. — Temperature as a function of distance from the AGB star in loga- 
rithmic scale used in our calculations for different wind acceleration radii. 
The temperature is ~ lO^K at lOAU where dust is formed. The crosses 
give dust t emperature estimates from IR observations of the symbiotic system 
Mira AB (Marengo et al.|2001^ . 




Fig. 6. — Schematic representation of the grid geometry in polar coordinates. 
The physical quantities are defined in the center of the cells. The system is 
centered on the primary and rotating in clockwise direction. 

loss from the disk around secondary. We used a grid with a 
fixed uniform size in the radial and azimuthal directions. The 
grid cells were sized to give square cells at the location of the 
secondary's orbit. To improve the resolution close to the sec- 
ondary additional levels of refinement were used to speed up 
the simulation. We introduce the mass of the secondary over 
1 orbital period to avoid the formation of strong shocks. The 
adaptive mesh refinement allows us to resolve a broad range 
of density contrasts close to the secondary. A diagram of the 
grid geometry is shown in Fig.|6] 

We assume that the wind is driven by radiation pressure and 
is composed of dust and a photoionized ideal monoatomic gas 
which are perfectly coupled. The wind acceleration mecha- 
nism is not considered in our model. The wind is accelerated 



at the dust condensation radius according to the Eq. 1 8 where 
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the full gravity of the AGB star is included. We include gas 
forces due to pressure and gravity from the AGB and the com- 
pact companion but the self-gravity of the gas is neglected. In 
our adiabatic runs, we assume that the gravity of the red giant 
is balanced by radiation pressure on the dust, so that it does 
not appear explicitly in our equations. Dust grains form in 
the AGB envelope, where pulsational shocks bring material 
from the atmosphere, and deposit it in the region of ~ 2 pho- 
tospheric radii. In this case the only force acting on the wind 
from the AGB star is the gas pressure, whereas in a binary 
the wind is additionally subject to the gravity of the compan- 
ion. Such an approximation leads to terminal velocities of 
the wind ranging between 10-40kms ', comparable to those 
obtained by Winters et al. (2000 1 for winds in which radiation 



pressure nearly balances the gravity of the wind source (see 
e.g. |Gawryszczak et al.^20 02). The effects of cooling are ac- 
counted for by a nearly isothermal equation of state in some 
of our runs, while in the adiabatic runs the polytropic index y 
is in the range 1 - 5/3. 

We neglect radiation pressure from the secondary. In the 
case of a white dwarf companion, for typical dust grain sizes 
the effect of radiation pressure on grains is negligible com- 
pared with gravity focusing. Explicit radiative heating or 
cooling terms are also neglected in our wind model. 

3.4.2. Equation of state 



The pressure is obtained by 



'RpT 



(19) 



where 'R is the gas constant, T is the temperature and ;U is the 
mean atomic weight. 

The temperature is approximately constant at the sec- 
ondary's position. Because of the gravitational field of the 
accretor, the disk aspect ratio calculated with respect to the 
/ib = H/\r - TbI can decrease considerably. This influences 
the flow in vicinity of the secondary and increases the amount 
of gas accumulated inside the Hill sphere. On the other hand 
we should expect a temperature increase in the vicinity of the 
secondary. For of a permanent disk hs can increase due to the 
compression of the gas flowing into the vicinity of the accre- 
tor. We use the following expression to calculate the sound 
speed 



GMa GMb 



((/za^-a)" + (fiBrB)") 



l/n 



(20) 



where ta = |r - TaI is the distance from the mass-losing pri- 
mary, tb - |r - TbI is the distance from the secondary, and a 
parameter n = 2 is used in our simulations. 

A disk with constant aspect ratio /ib can be formed around 
the mass gainer with a Keplerian angular velocity 



GMn 



(21) 



3.4.3. Model parameters 

The main parameters of the models are the wind tempera- 
ture that determines its terminal velocity, the orbital parame- 
ters (binary separation and orbital period for a circular orbit), 
mass loss rate from the AGB star and the mass ratio defined 
by 



TABLE 1 

Model parameters of a long-period and close symbiotic binaries. 



Parameter long-period 

Mass of primary Ma (g) 2.387 X 10^'' 

Mass of secondary Mg (g) 1.193 X 10^^ 

Radius of primary «a (cm) 2.394 X lO'" 

Radius of secondary (cm) 1.496 X lO''' 
Wind temperature T^wind 

(K) 1000 

Binary separation d (cm) 1.047 X lO''' 
Period Pb (yr) 436.5 



short-period 
2.387 X lO" 
1.193 X 1033 
2.394 X 10'" 
1.496 X 10'" 

1200 
2.394 X 10'" 

40.9 



where Ma is the mass of the primary mass-losing star and Mb 
is the mass of the wind-accreting compact star 

Here we describe the model parameters of the runs dis- 
cussed in Section |4] The indexing of the models described 
below increases with increasing dust condensation radius to 
facilitate reference to models. The average wind terminal ve- 
locity from an AGB star is around 15 kms within a range 
of 10 kms In the following, we produce models having a 
terminal velocity for a single star between 10 and 25 kms"'. 
The binary separations d vary between 10 and 70 AU. Our 
base resolution with 3 levels of refinement corresponds to a 
cell size of ~ 0.004 close to the accretor Each simulation 
with our standard resolution requires about 30 hours of CPU 
time on Pentium 4 processors. 

The inner Lagrangian point in model M-4 is at 40 AU from 
the mass-lossing star, which is close to the inner boundary of 
the domain. In this case we expect to see some of the effects 
of RLOF in addition to that of wind accretion. 

The parameters in these simulations are closely related to 
two symbiotic systems Mira AB and R Aqr (see Table[T]). The 
wide symbiotic system Mira AB (o Ceti) consists of a pulsat- 
ing AGB primary (Mira) and a compact secondary at an Hip- 
parcus distance of ~ 128pc. Mira stars are long-period vari- 
able stars which evolve along the AGB and present large stel- 
lar pulsations. The primary's radius is found to be ~ 270Rq. 
The nature of the companion is uncertain and has been sug- 
gested to be a white dwarf or a late-type, magnetically active 
main-sequence dwarf. HST observations of Mira AB have 
shown that it is a strongly interacting system with mass trans- 
fer from the primary towards the compact companion. Its op- 
tical through UV spectrum appears to be dominated by emis- 
sion from an accretion disk that presumably is accumulated 
from Mira A's massi ve wind (Rei mers & Cassat ella 1985] 
Karovska et al.|I997) l. Recent observations of symbiotic bi- 
naries in wavebands ranging from the IR to the X-rays sug- 
gest that wind focusing can form a s tream flow even in wide 
systems ( Karovska et al.|[T997 2005). The system has large 
variations in luminosity with well-defined pulsation periods. 
X-rays from the Mira AB system arise either from magneto- 
spheric accretion of wind material from Mira A onto Mira B, 
or from coronal activity associated with Mira B itself, as a 
consequence of accretion-driven spin-up. 

R Aqr is a close symbiotic binary system (at a distance of 
~ 200pc) composed of a mass-losing Mira variable and a 



white dwarf surrounded by an extended nebula ( Solf & Ul- 
rich 1985| l. The orbital period is not well determined with 



values in the literature ranging between 2 and 44 yr d Willson] 
etal.|1981||Hinkle eta l. 1989 ; Mcintosh & Rustan 2007ji. The 
companion is surrounded by an accretion disk which gives 
rise to a jet with bipolar outflow. The primary almost fills its 
Roche lobe during periastron, although mass flow towards the 
companion is increased due to wind focusing during the rest 
of the orbital motion ( |Gromadzki & Mikolajewska|2008| l. 
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4. RESULTS AND DISCUSSION 

The wind parameters and mass ratio of our models of wind 
accreting symbiotic binaries are given in Table [2] We will 
consider mainly the simulations with a locally isothermal 
equation of state in the following discussion. The adiabatic 
runs will be considered in Sec. 14.2.21 

4. 1 . Roche lobe overflow 

First, we study the accretion in the RLOF case to check our 
code and compare the results with the wind accreting cases 
described below in Sec. 4.2 Also, in some symbiotics RLOF 



is likely to occur during periastron if their orbits are highly 
elliptical, for e xample in R Aqr where the stellar separation 
becomes 8 AU ( [Gromadzki & Mikolajewska|2008 ). The gas 
from the primary star is assumed to flow through the inner 
equilibrium point with an initial velocity between 0.01 - 0.1 
of the orbital velocity. We run our simulations for 5 orbital 
periods, when the system has reached a quasi-steady state. 

Our simulations show that most of the flow is efficiently 
accreted onto the star as it turns around and forms a stable 
disk with counterclockwise rotation and elliptic streamlines. 
The gas is accelerated by the companion's gravity and Cori- 
olis force and becomes supersonic. The Li stream cannot hit 
the companion star due to the Coriolis force and penetrates 
into the disk forming a hot spot that heats up the outer disk in 
our adiabatic runs with y = 5/3. 

In Fig. |7] we show the density contours for a system with 
separation of 70AU and Li Lagrangian point located at ~ 
30AU from the accretor using two diff'erent equations of state. 
Mass is removed inside 0.5 Hill radius in our simulations to 
achieve a quasi-stable state after several orbital periods. The 
angle of deflection of the stream roughly agrees with the ana- 
lytical estimates of |Lubow & S hu ( 1975 ). Fig.[8]shows the ve- 
locity vectors in the inertial frame for one of the RLOF runs. 

The flow patterns are non-axisymmetric due to the shock 
waves that form when gas rotating in the disk collide with 
the stream. This kind of patterns have been observed using 
Doppler tomography in the disk around IP Pegasi (Steeghs 

2000| l. The angular momentum 



et al. 1997 Makita et al. 



transfer in the disk is more efficient due to the presence of 
shocks and accretion may take place without turbulent vis- 
cosity. There is some mass loss in our grid in the case of 
reflective boundaries at the accretor's surface. 

Since cooling effects due to radiation are not included in 
our equation of state we use a locally isothermal model where 
the temperature decreases as r ' from the center of the sec- 
ondary. In this case the stream shocks are stronger than in the 
adiabatic runs (right panel in Fig.[7]i. 

4.2. Wind accretion 

Our results indicate that wind accretion in symbiotic bina- 
ries can result in RLOF-like accretion via a stream flow. Sta- 
ble accretion disks can form f or winds with ma ss loss ranging 
between 10 - 10 ^'Mgyr ' ( |Ryde et al.|2000| l and distances 
between 16 - 70AU. 

The mass accretion rates onto the secondary are of the or- 
der of the predicted ones from Bondi-Hoyle theory when the 
inner boundary of our grid is less than 20AU. Here, we will 
focus on slow winds with a relative speed at the position of 
the accretor of 6 - 25km s"'. Mass accretion ratios for differ- 
ent values of the accretion radius r^cc inside the Hill radius in 



Eq. 15 give very similar results after we reach a stable con- 



TABLE 2 

Models with separation 70 AU and mass ratio ( 

FORMATION RADIUS. 



: 2 FOR DIFFERENT DUST 



Model 


Dust formation radius 
(AU) 


Wind temperature 
(K) 


q 


hB 


M-1 


10 


900 


2 


0.1 


M-2 


20 


490 


2 


0.1 


M-3 


30 


340 


2 


0.1 


M-4 


40 


300 


2 


0.1 


M-5 


30 


300 


2 


0.05 



around the accretor, the accretion ratio is consistent with the 
lower resolution simulations. 

For a fast wind we expect an accretion rate of th e order 
of the Bondi-Hoyle accretion ( Bondi & Hoyle|1944[ l. In this 
case we find a quasi-spherical wind morphology. For wind 
acceleration radii between 30 to 40AU the accretion rate is a 
few times the Bondi-Hoyle rate with variations with an ampli- 
tude of about 10%. The wind presents a broad range of den- 
sity contrasts in our computational domain. Bow shocks are 
formed close to the secondary in these models at late times. 
It has been suggested that UV radiation in symbiotic systems 
is originated in those shocks ( Theuns & Jorissen 1993). In 
Table|2]we show the parameters of our models that reproduce 
closely the binary system Mira AB for different wind accel- 
eration radii. For typical mass loss rates from an AGB star 
( Ryde et al.|2000| l we compute an accretion rate on the com- 
panion of ~ 10" - 10"'^ M© yr"' via stream flow accretion. 
The X-ray luminosity of the accretor can be estimated using 



130 



(23) 



which corresponds to ~ 3 x 10^" - 3 x 10^^ ergs" 
ing a radius - 0.02Ro for the compact object in our M-3 
model. 



In Fig. 



12 we compare the accretion ratio in our mod- 
els with the Bondi-Hoyle accretion given by Eq. [3] Table |3] 
shows accretion rates onto the compact object divided by the 
mass loss rate of the AGB star as a function of time. Accre- 
tion ratios become quasi-stable after half of an orbit of the 
system. However, in the Bondi-Hoyle formula the velocity 
is the upstream wind velocity at infinity, while in our binary 
simulations we use the relative wind velocity at the secondary 
position. Moreover, the plane parallel flow assumption in the 
Bondi-Hoyle theory is not valid for our wind model at close 
distances. 

A small change in the wind parameters can lead to changes 
of a few times in the accretion rate when an accretion disk is 
formed around the secondary. In particular, it could explain 



the increase of extinction in symbiotic binaries ( Gromadzki & 
Mikolaj ewska|2008 1 as a change in the accretion flow. Those 



variations can occur in a fraction of the orbital period and may 
also explain Nova outbursts in objects such as RS Ophiuchi 
(IWalder et al.T2008]l. 



In Fig. 13 the azimuthally averaged density distribution 
close to the accreting star for model M-3 is shown. The den- 
sity decreases exponentially up to 5 AU in a slightly eccen- 
tric disk and is likely to be optically thick. This is consistent 
with the size of the disk around Mira B estimated by [Ireland] 



et al. ( 2007| l. We show the flowlines after 1 orbital period 



figuration. In our runs with 4 additional levels of refinement 



in Fig. 14 The flowlines are calculated integrating the equa- 
tion of motion from the line joining the stars using positive 
and negative time-steps. When the wind speed is of the order 
of the escape velocity from the AGB star the flow is focused 
by the gravitational pull from the secondary and the accretion 
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Fig. 7. — Density contours in tlie orbital plane of an RLOF simulation for an adiabatic run with y : 
has a separation of 70AU and a mass ratio q = 2. 



1.1 (left) and a locally isothermal run (right). The system 
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Fig. 8. — Density contours in the orbital plane in polar coordinates for the 
RLOF simulations after 1 orbit. Velocity vectors are shown in the inertial 
frame . 

rate is enhanced by a factor of a few with respect to standard 
Bondi-Hoyle accretion. 

We show the radial velocity map for a model with wind 
temperature 10^ K and inner boundary at 20 AU with over- 
plotted grid structure in Fig. 15 The refinement criterion used 
by PARAMESH is based on the second derivative normalized 
by the gradient over one cell (Lohner 1987). The blocks are 
approximately square close to the secondary and contain 8x8 
computational cells. 

The orbital speed in the disk is close to Keplerian with a 
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Fig. 9. — Temperature distribution in logarithmic scale after 6 orbits for a 
binary system with separation of 70AU. 



Steeper radial profile. In Fig. 16 we plot the averaged az- 



imuthal velocity around the secondary. 

4.2.1. Influence of the wind temperature 

We tried several models to check the dependence of the 
wind parameters on the formation of a stream flow between 
the stars. The dust is believed to be condensed at a temper- 
ature of ~ 1500K which corresponds to a distance of 8 AU 
from the center of the primary (Ivezic & Elitzur|1997)l. When 
the wind acceleration radius is of the order 50% of the Hill 
radius a density enhancement is formed focusing 10% of the 
slow wind inside the Roche lobe. The sound speeds in our 
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Fig. 10. — Density contours in logaritlimic scale for a separation of 70AU at several times. From left to right the inner boundary where the wind is accelerated 
is 10,20,30 and 40AU. The model parameters are presented in Table[2] In the outermost right column the envelope of the primary is close to the Hill surface. 
Three snapshots are shown for times / = 0.25,0.5 and 1 orbital periods. The temperature of the wind decreases roughly as r"' where r is the distance from the 
primary. 

models are in the range 3 - 10km s ' which are of the same 
order of the wind speed relative to the accretor When the 
average velocity of the wind at the inner Lagrangian point is 
larger than the orbital velocity of the system, the stream does 
not form and the accretion ratio is close to Bondi-Hoyle esti- 
mate. 



In Fig. 10 we show the density contours for a binary sepa- 
ration of 70AU for several values of the dust formation radius 
where the wind is accelerated. In these simulations we include 
the full gravity of the primary. 

For a dust formation radius of 5AU the wind geometry is 
almost spherical. A bow shock forms as the wind is decel- 
erated by the gravity from the companion and the accretion 



ratio is close to the Bondi-Hoyle value. A small accretion 
disk is formed around the secondary. As we move the dust 
formation radius outwards the temperature for dust conden- 
sation decreases according to Fig [5] When we place the dust 
formation radius at lOAU the flow structure is clearly elon- 
gated and the accretion rates are a few times larger than the 
Bondi-Hoyle theoretical estimate. This is consistent with an 
agglomeration of dust particles at several stellar radii above 
the photosphere in AGB stars in order to have focused wind 
RLOF-like mass transfer via stream flow in wi de symbiotic 
systems such as Mira AB (Karovska et a l.|2()05 i 



For the cases where the dust formation is < 20AU the wind 
velocity is lower than the escape velocity in our model. The 
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Fig. 1 1 . — Density maps and velocity vectors in the orbital plane in loga- 
rithmic scale for a separation of 70AU. Three snapshots are shown for times 
t = 0.5, 1 and 1.5 periods after the secondary is introduced. 
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Fig. 12. — Accretion ratio, defined as the accretion rate divided by the 
mass-loss rate from the AGB star, as a function of time for our model M-5 
(see Table[2j. 



TABLE 3 

Accretion ratios as a function of time for the simulation shown in Fig. 1121 



Time (orbits) -P^l 



0.1 


0.0187499146615 


0.2 


0.032895734525 


0.3 


0.0405163717933 


0.4 


0.0457114513595 


0.5 


0.0486647436321 


0.6 


0.052120162122 


0.7 


0.0518595175058 


0.8 


0.0545558407777 


0.9 


0.0543526731365 


1.0 


0.0550009176397 


1.1 


0.0557649810422 


1.2 


0.0573139141097 


1.3 


0.0579923479442 


1.4 


0.0587394563781 


1.5 


0.0588350261856 


1.6 


0.058462177135 


1.7 


0.059570725981 


1.8 


0.0612462210888 


1.9 


0.0612788117294 


2.0 


0.0605904603547 



wind opposite from the companion is not able to overcame 
the gravity and forms a bow shock. This does not affect the 
formation of a tidal stream and our calculation of the mass 
accretion rate onto the secondary. The escape velocity from 
the primary is ~ 10km s"' at the dust formation location in our 
model M-2 (second column in Fig.fTOl). In the slow wind case, 
numerical simulations show intricate flow patterns, therefore 
it is not possible to apply directly the results from Bondi- 
Hoyle accretion theory. 

In our model M-5, we have changed the temperature profile 
around the secondary The accretion ratios averaged over one 
orbital period agree within 5% in our models M-3 and M-5. 

4.2.2. Modifying y 

We have varied the poly tropic index between y - 1 - 5/3. 
The density distributions are similar within 10% for the sim- 
ulations with q ^ 2, d ^ IQAU and y = 1,1.2,5/3. The 
variation of the polytropic index intends to bracket more real- 
istic cases that include radiative losses in the equation of state. 



Numerical simulations of wind accretion 



11 



t= 1.000 



-4.5 



-5.0 - 



CO 

o 



-5.5 - 



-6.0 




i-[AU] 

Fig. 13. — Averaged radial density profile of the disk around the accreting 
secondary in logarit hmic scale in computational units after 1 orbit for the 
model shown in Fig. 1 1 1 1 The distance in the horizontal axis is given in AU 
from the center of the secondary. 
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Fig. 14. — Streamlines in the corotating frame after 1 orbit overplotted on 
the dens ity contours in logarithmic scale for the same model as shown in in 
Fig. |ll| The flowlines around the secondaiy are eccentric. 

In the adiabatic simulations the accretion ratio varies between 
8 - 10% and is several times larger than in the Bondi-Hoyle 
scenario. 

4.2.3. Variation of the orbital parameters 

The orbital velocity in our models varies between 5 - 
10km s ' while the wind velocity at the surface of the AGB. 
is of the order 2 - 10km s"' . 

the density contours of a system with separation 



In Fig. 17 



of 16AU and mass-losing surface located at 8AU from the ori- 
gin is shown as a function of time. The central star has a mass 
of 1.2Mq and an accreting companion of 0.6Mq orbits on a 
fixed cuxular orbit. An accretion disk is observed around the 
secondary with Keplerian velocity profile and a narrow stream 




Fig. 15. — Radial velocity contours in the corotating frame after 0.5 or- 
bits with grid structure using 3 levels of refinement. Each block has 8x8 
computational zones. 




Fig. 16. — Averaged orbital velocity profile of the accretion disk around the 
secondary for the model shown in Fig. jl 1| The distance in the horizontal axis 
is given in AU from the center of the secondary. 

forms between the stars. The accretion flow focused onto the 
secondary is about 15% of the unperturbed wind mass loss 
which is several times more efficient than the models for bi- 
nary separation of 70AU. Mass is removed from the accretion 
disk after each time-step to avoid the accumulation of mass 
and formation of strong shocks according to Eq. 15 (Gunther 
& Kley|2002| . 



4.2.4. Convergence tests 

We have rerun simulation M-3 at a twice the linear reso- 
lution of our reference simulations with three levels of re- 
finement, and performed an additional run at the same base 
resolution with four additional levels of refinement. In our 
high resolution runs, the gravitational softening e is half the 
size of the gravitational softening in the low resolution runs. 
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Fig. 17. — Density contours in the orbital plane with overplotted velocity 
vectors in the corotating frame for a system with separation of 16AU after 
0.25, 0.5 and 1 orbital periods. 



The density distributions in the higher resolution simulations 
agree with our previous results within 5% with larger oscilla- 
tions. Therefore, we conclude that the resolution of our simu- 
lations is acceptable and does not affect our results. 

5. CONCLUSIONS 

We have constructed a hydrodynamical time-dependent 
model of a wind from an AGB star which is gravitationally 
focused by a close companion. The model has several free 
parameters. The mass ratio q of the system, the binary sep- 
aration d and the eccentricity e characterize the binary con- 
figuration. The effective temperature on the surface of the 
mass losing star, the radius of the AGB star where the dust is 
formed and the mass loss describe the properties of the wind 
from the mass-losing star We have improved on previous hy- 
drodynamic simulations by using an equation of state where 
the temperature distribution depends on the distance from the 
primary star and the accreting companion. 

For each binary model we integrate the equations of hydro- 
dynamics until the density distribution of the gravitationally 
focused wind reaches a stationary state. The general effect 
of the secondary is to deflect the AGB wind toward the or- 
bital plane, and as a result the relaxed density distributions 
present accretion disks of various characteristics. Accretion 
via a gravitationally focused wind has implications for a broad 
range of symbiotic and wide binary systems and can explain 
the formation of chemically peculiar stars. A tidal stream is 
observed in our runs for binary separations up to ~ IQAU 
for slow winds and nearly isothermal flows. Wind focusing 
weakens with increasing orbital separation because of the in- 
creasing ratio of the wind velocity at the secondary's position 
to the orbital velocity of the system. 

All our models develop a high density region in the form 
of a Keplerian thin disk around the secondary with a mass of 
about 10""* Mq and an exponentially decreasing density pro- 
file. The temperature profile around the accretor was fixed 
in our simulations and corresponds to a disk scale height of 
~ 0.05. Winds from evolved stars in binary systems show 
complex flow patterns and are far from being spherically sym- 
metric. A thin accretion disk around Mira B has been inferred 



from spectroscopic UV observations (e.g. Reimers & Cas- 
[satella 1985; Karovska et al. 1997 ) and interferometric imag- 
ing (Ir eland et al.„2007^ . The formation of a stream flow is 
strongly dependent on the wind temperature at the dust for- 
mation radius which determines the wind velocity at the inner 
Lagrangian point. For the winds parameters from a AGB star 
(Mira type) we obtain stable disks with masses of the order 
10""* of the stellar mass for a range of orbital separations. The 
accretion rate onto the secondary is almost constant once the 
accretion disk is formed and we reach a quasi-stable state. We 
conclude that the main parameter that determines the forma- 
tion of a stream and enhanced accretion rates compared to 
Bondi-Hoyle's rates is the wind velocity at the dust forma- 
tion radius. Nevertheless, the Bondi-Hoyle estimate does not 
include pressure forces and is not directly applicable to our 
models of wind accretion in symbiotic binaries. 

The presence of a long-lived disk and stream flow is sensi- 
tive to the wind parameters, in particular of the dust accelera- 
tion radius where we place the inner boundary of our grid. In 
case that the wind parameters change with time the accretion 
rates can be expected to change in relatively short time-scales 
due to rapid changes in the flow patterns and the formation 
of a stream flow. This can explain the variability in accre- 
tion rate and transitions from a quiescent to an active state in 
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binary systems ( Bisikalo et al.|[2002 l. As the secondary in- 
creases in mass there are several possible outcomes for the 
system. The angular momentum loss in symbiotic binaries 
may be large enough to imply a shrinking of the orbit and a 
possible merger that can l ead to the formation of B arium stars 
or Type la supernovae ( ^Theuns et al. 1996 Jahanara et al., 
|2005] l. 

The massive wind of an AGB star can be considerably de- 
flected towards the orbital plane by the gravitational interac- 
tion with a companion star. The X-ray luminosity powered by 
accretion onto the mass gainer can photoionize the wind and 
reduce its ability to be accelerated by radiation pressure. The 
lower radiative driving can result in a reduced wind velocity 
near the companion, which further increases the the accretion 
efliciency of the companion. In the case of rapid accretion 
rate changes a symbiotic binary may undergo outbursts such 
as those observed in Mira AB ( |Karovska et al.|2005] l. A de- 
tailed comparison of the results from the numerical models 
with multi-wavelength observations of Mira AB will be pre- 
sented elsewhere. 

We have adopted a simple a simple prescription for the 
equation of state and wind acceleration mechanism. In the 
future our calculations can be improved by using a more de- 
tailed energy balance for the gas going beyond the isother- 



mal approximation including radiation cooling and heating 
due to shock formation. Radiative feedback is also expected 
to be important in the context of symbiotic systems and X- 
ray binaries ( [Edgar & Clarke 2004|. We plan to perform 
3-dimensional simulations to compare our numerical results 
with high resolution X-ray observations of the Mira AB sys- 
tem. 
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